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1.  Introduction 


Future  military  engagements  will  require  weapons  systems  exhibiting  longer  range  and  greater 
accuracy.  One  of  the  technologies  under  investigation  to  achieve  these  goals  is  the 
electrothermal-chemical  (ETC)  propulsion  concept,  shown  schematically  in  figure  1.  In  the  ETC 
gun,  energy  which  is  stored  either  in  batteries  or  in  a  rotating  device  is  converted  on  demand  into 
an  electrically  generated  plasma  (resulting  from  the  ablation  of  polyethylene  material  in  a 
capillary),  which  is  then  injected  into  the  chamber  in  a  howitzer  or  gun.  This  plasma  energy  is 
used  to  ignite  the  chemical  propulsion  charge  (e.g.,  solid  propellant)  as  well  as  enhance  gun 
performance  by  taking  advantage  of  a  number  of  unique  plasma  characteristics.  For  example,  a 
low-density  plasma  jet  can  more  efficiently  ignite  charges  of  high  loading  density  (HLD),  can 
control  propellant  mass-generation  rates  (7),  can  reduce  propellant  charge  temperature  sensitivity 
(i.e.,  the  variation  of  gun  performance  with  changing  ambient  temperature  [2,  3])  and  can  shorten 
ignition  delay  (i.e.,  the  time  interval  between  firing  of  the  igniter  and  ignition  of  the  propellant 
[7]).  Plasma  igniters  also  eliminate  the  conventional  chemical  igniter  and  can  thus  enhance  the 
safety  aspects  of  the  overall  gun  propulsion  system.  All  of  these  observations  have  a  significant 
effect  on  the  ballistics  of  ETC  gun  systems  and  can  lead  to  a  useful  improvement  in  gun 
performance  and  accuracy. 


Figure  1.  ETC  gun  concept. 


Research  has  been  carried  out  on  a  number  of  aspects  of  plasma  and  plasma/propellant 
interactions  (5,  6).  Because  the  plasma  is  at  a  temperature  (typically  >10,000  K)  that  is 
considerably  higher  than  conventional  chemical  igniters  (3000  K),  the  radiation  properties  of  the 
plasma  have  also  been  considered.  The  high  plasma  temperature  leads  to  radiation  effects  nearly 
100  times  greater  than  that  of  chemical  igniters  (i.e.,  a  T4  effect)  (5).  Such  radiation  could  lead 
to  significantly  different  temperature  profiles  within  the  propellant,  causing  changes  in 
combustion  rates.  Another  significant  characteristic  of  the  plasma  is  a  much  lower  density  than 
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the  gases  generated  by  a  chemical  igniter.  The  lower  density  plasma  impacts  the  convective 
heat  transfer  to  the  propellant,  as  the  plasma  moves  through  the  grains,  as  well  as  the  velocity 
and  mode  of  flamespreading  within  a  propellant  bed.  It  has  been  suggested  that  energy  transport 
by  convection  may  be  as  important  as  radiation  transport  in  plasma-propellant  interaction  (PPI) 
(7,  5). 

The  effects  previously  described  can  lead  to  significant  changes  in  ballistic  behavior  and  useful 
improvements  in  gun  performance  once  an  understanding  of  the  underlying  physical  mechanisms 
is  achieved.  To  this  end,  the  U.S.  Army  Research  Laboratory  (ARL)  has  begun  a  comprehensive 
study  on  the  interaction  of  the  plasma  efflux  from  a  ETC  igniter  with  solid  propellant  grains. 

The  goal  of  this  work  is  to  elucidate  the  relevant  physical,  mechanical,  and  chemical  mechanisms 
that  underlie  the  observed  ballistic  effects.  Various  aspects  of  the  experimental  program  for  PPI 
are  described  elsewhere  (9).  The  first  phase  of  the  modeling  effort  in  support  of  this  PPI  project 
involves  a  time-accurate  computational  fluid  dynamics  (CFD)  code  that  has  been  written  to 
include  high-temperature  thennodynamics,  variable  specific  heats  and  transport  properties 
(viscosity  and  thermal  conductivity),  and  finite-rate  (nonequilibrium)  chemical  kinetics  (the 
mechanism  is  described  in  reference  [10]).  A  separate  capillary  model  described  by  McQuaid 
and  Nusca  {11,  12)  supplies  boundary  conditions  for  the  CFD  code  in  terms  of  the  physical  and 
chemical  properties  of  the  plasma  capillary  efflux.  Validation  of  the  capillary  model  and  the 
CFD  code,  including  coupled  chemistry,  uses  a  series  of  experiments  wherein  a  plasma  jet  is 
generated  from  a  plasma  cartridge;  pressures  in  the  resulting  unsteady  flowfield  are  measured 
using  a  probe. 


2.  Description  of  Open-Air  Plasma  Jet  Experiment 


References  (5)  and  {13)  provide  a  full  description  of  the  experiment  with  a  summary  given 
herein.  The  general  view  of  the  plasma  generator  is  shown  in  figure  2.  The  plasma  capillary 
consists  of  standard  polyethylene  liner  with  dimensions  of  102  mm  in  length  and  6.35  mm  in 
diameter.  The  nozzle  consists  of  a  circular  tube  of  2.9  cm  in  length  and  1 .3  cm  in  diameter  that 
is  directly  adjacent  to  the  end  of  the  plasma  capillary  followed  by  another  circular  tube  of  1 .9  cm 
in  length  and  3.2  cm  in  diameter.  As  a  result,  the  plasma  efflux,  generated  by  the  capillary,  is 
expanded  through  a  step-configuration  nozzle  before  emerging  into  the  open  air.  A  series  of 
probes  is  positioned  in  the  path  of  the  plasma  efflux;  typically,  a  single-pressure  probe  is  used 
(figure  2),  but  ongoing  tests  involve  heat  flux  probes  as  well  {13).  To  protect  these  probes,  the 
plasma  capillary  is  mounted  behind  a  steel  plate.  The  surface  of  this  plate  forms  the  upstream 
boundary  for  the  numerical  simulation,  with  a  hole  in  the  plate  fonning  the  inlet  to  the 
computational  domain. 

The  pressure  probe  is  placed  1 5  cm  away  from  the  surface  of  the  plate  and  aligned  with  the 
centerline  of  the  cartridge.  The  data  acquisition  system  was  chosen  for  maximum  frequency 
response  (>1  MHz).  The  pressure  gage  was  rated  at  1.5  MHz  and  0.2  ps.  Care  was  exercised 
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PFN 


Figure  2.  Schematic  of  plasma  generator  and  pressure  probe  also  showing  computational  domain. 

in  shielding  the  gage  from  direct  interaction  with  the  plasma  as  were  the  cables  and  amplifiers 
(strong  evidence  of  charged  particle  impact  on  the  unshielded  gage  was  found  [5]).  Even  so, 
post-firing  examination  of  the  pressure  gage  indicated  significant  amounts  of  particulate  in  the 
plasma  effluent  (the  plasma  is  generated  by  ablation)  that  contributed  to  noticeable  variability  in 
pressure  histories  from  shot  to  shot. 

The  electrical  power  supplied  to  the  plasma  generator  used  a  five-stage,  pulse-forming  network 
(PFN)  that  could  store  up  to  130  kJ  (see  reference  [5]).  The  five  stages  could  be  charged  to 
different  voltages  and  fired  at  different  times  to  obtain  the  desired  power-pulse  shape.  The 
electrical  data  were  recorded  using  a  Rogowski  coil  for  load  current  and  a  voltage  divider  and 
current  transformer  for  load  voltage.  These  data  were  used  to  generate  load  current,  voltage, 
power,  energy,  and  impedance  as  a  function  of  time.  Figure  3  shows  the  power  and  amplitude 
variation  for  a  typical  experiment.  Given  the  amplitude/time  and  the  physical  characteristics  of 
the  capillary,  the  capillary  model  (11, 12)  generates  the  range  of  density,  velocity,  pressure, 
temperature,  and  species  distributions  shown  in  figures  4-6. 

Note  that  figure  3  shows  a  pulse  forming  network  (PFN)  pulse  width  of  -0.45  ms,  but 
figures  4-6  show  efflux  histories  spanning  -0.275  ms.  The  current/time  data  up  to  0.027  ms 
and  between  0.40  and  0.45  ms  were  ignored  by  the  capillary  model;  the  relatively  small  current 
levels  during  these  times  caused  a  divergence  in  the  capillary  model  because  a  relatively  small 
amount  of  plasma  was  being  generated.  As  a  result,  the  efflux  histories  were  shifted  to 
zero-time,  and  the  CFD  computation  was  started  with  efflux  properties  that  characterized  an 
abrupt  rather  than  a  more  gradual  start  to  the  efflux;  the  computed  time-of-arrival  of  the  jet  on 
the  pressure  probe  will  be  shorter  than  measured  (see  section  5).  A  more  amenable  resolution  to 
this  problem  is  being  formulated  for  future  simulations.  Figure  6  shows  that  there  are  five 
prominent  species  in  the  plasma  efflux  among  a  total  of  3 1  species  that  were  considered  in  the 
capillary  model  (11, 12)  and  38  species  in  the  CFD  simulations  (7  species  were  added  for 
reaction  with  air).  For  a  complete  list  of  species,  see  table  1. 
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Figure  4.  Efflux  density  and  velocity  histories. 
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Figure  5.  Efflux  pressure  and  temperature  histories. 
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Figure  6.  Efflux  species  histories. 


hed  Line) 

15500 

-  15000 

■  14500  ^ 

■  14000 

-  13500 

-  13000 

■  12500 
1.3 


TEMPERATURE 


Table  1.  Species  concentrations  (moles/m3)  at  0.22  ms  and  four  locations  along  the  centerline. 


Species 

Electrons 

C 

c+ 

C++ 


ho2 

h2o2 

hno2 

no2 

co2 

o3 

Density  (mol/m3) 
Pressure  (atm) 
Temp  (K) 
‘Quantity  is  <1.0  x  10 


0.0022 

0.0061 

0.0017 

0a 

1.07E-06 
3.32E-06 
1.90E-06 
9.48E-1 1 
9.16E-12 
2.58E-10 
3.27E-11 
2.12E-06 
9.12E-07 
0.0152 
5.40E-04 
2.40E-06 
4.28E-06 
5.43E-07 
9.25E-08 
5.34E-09 
4.98E-11 
1.14E-1 1 


3.31E-16 

0a 

1.33E-07 

4.32E-09 

6.80E-11 

1.33E-11 

0a 

2.72E-17 


0a 

0.0274 

31.7 


X  =  1  cm 


0.0026 

0.0048 

0.0019 

0a 

1.09E-06 
2.38E-06 
1.99E-06 
5.29E-11 
8.63E-12 
3.41E-10 
2.59E-06 
1.27E-06 
8.16E-07 
0.0126 
7.61E-04 
2.78E-06 
3.52E-06 
7.46E-07 
7.60E-08 
7.62E-09 
3.78E-11 
1.44E-1 1 
1.03E-16 
2.68E-11 
1.95E-18 
2.70E-06 
6.24E-09 
1.10E-10 
1.67E-1 1 
4.98E-16 
8.96E-13 
1.08E-15 


2.72E-21 

0a 

0.0194 

21.1 


X  =  9  cm 


0.0021 

0.0038 

0.0015 

0a 

8.53E-07 
1.19E-06 
1.55E-06 
4.32E-1 1 
6.79E-12 
5.42E-10 
2.07E-06 
1.03E-06 
6.45E-07 
0.0099 
5.80E-04 
2.15E-06 
2.82E-06 
5.71E-07 
6.03E-08 
5.81E-09 
3.02E-1 1 
1.11E-11 
5.80E-17 
2.14E-1 1 
3.82E-12 
2.15E-06 
5.02E-09 
1.09E-10 
1.28E-1 1 
6.26E-10 
7.39E-10 
1.06E-15 


1.21E-20 

0a 

0.0153 

22.7 


X  =  15  cm 
(Probe  Tip) 

0.0435 

0a 

0.0327 

0a 

2.38E-07 
6.66E-07 
4.27E-07 
1.41E-1 1 
1.99E-12 
0.0017 


4.05E-07 

1.95E-07 

0.0462 

0.0429 

5.57E-07 

9.22E-05 

1.34E-07 

7.17E-08 

1.33E-09 

4.45E-12 

2.72E-12 

3.92E-12 

5.17E-08 

0.0061 

0.0767 

0.0014 

2.02E-04 

3.17E-12 


1.46E-08 


1.35E-08 

0a 

0.0894 

62.8 
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3.  Previous  Modeling  Efforts 


ARL,  in  cooperation  with  the  Naval  Research  Laboratory  (NRL),  had  conducted  a  computational 
study  of  open-air  plasma  discharges  from  1997  to  1998  (14-16).  This  computational  study  was 
based  on  the  NRL  CFD  code  FAST3D  (1 7)  and  sponsored  by  the  Department  of  Defense 
Common  High-Performance  Computing  Software  Support  Initiative  program.  The  ARL  open- 
air  plasma  firing  was  chosen  as  a  test  case  for  the  FAST3D  CFD  code  because  the  plasma  jet  is 
of  a  very  high  velocity  and  temperature,  but  of  a  very  low  density.  The  plasma  jet  was  assumed 
to  be  an  ideal  neutral  single-component  gas  with  density,  temperature,  ratio  of  specific  heats  (y), 
and  velocity  characteristics  identical  to  the  plasma  discharge.  The  multicomponent  nature  of  the 
plasma,  consisting  of  various  neutral  species,  ions,  electrons,  and  its  chemical  reactivity  was  not 
considered  in  those  studies.  The  augmentation  of  plasma-air  mixing,  caused  by  viscosity  and 
turbulence,  was  also  excluded  in  these  studies.  The  time-dependent  gas  properties,  along  with  jet 
velocity  and  y,  were  determined  using  a  plasma  capillary  code  developed  at  ARL  by  Powell  and 
Zielinski  (18).  Overall,  these  CFD  simulations  captured  the  major  features  of  the  plasma  jet  as 
photographed  by  White  et  al.  (5).  Alternating  bright-dark  structures  in  the  photos  were  identified 
as  the  shock  structures  in  the  jet.  A  precursor  shock  caused  by  the  high-velocity/low-density 
plasma  effluent  (a  feature  not  easily  observed  or  photographed)  was  captured.  Comparisons 
between  computed  and  measured  pressures  in  the  jet  were  encouraging  but  pointed  to  the  need 
for  a  variable  y  multicomponent-reacting  flow  model. 


4.  Multispecies  Reacting  Flow  CFD  Code 


The  high-temperature,  nonideal  chemically  reacting  gas  flowfield  within  the  capillary  efflux  jet 
is  numerically  simulated  using  CFD.  The  NSRG2  code,  written  by  Nusca  (19-21),  solves  the 
two-dimensional  (2-D)/axisymmetric,  unsteady  real-gas  Navier-Stokes  equations,  including 
submodels  that  represent  finite -rate  (nonequilibrium)  chemical  reactions,  multispecies  diffusion, 
as  well  as  variable  specific  heats,  viscosity,  and  thermal  conductivity.  These  partial  differential 
equations  are  cast  in  conservation  form  and  converted  to  algebraic  equations  using  a  finite- 
volume  formulation.  Solution  takes  place  on  a  mesh  of  nodes  distributed  in  a  zonal  fashion 
throughout  the  flowfield  such  that  sharp  geometric  corners  and  other  details  are  accurately 
represented.  The  conservation  law  form  of  the  equations  ensures  that  the  end  states  of  regions  of 
discontinuity  (e.g.,  shocks)  are  physically  correct  even  when  smeared  over  a  few  computational 
cells.  The  Navier-Stokes  equations  for  2-D/axisymmetric  (x,y  coordinates),  reacting  (N  species), 
and  unsteady  (time,  t)  flow  are  written  in  nondimensional  form  (19).  The  dependent  variables 
are  density  (p),  velocity  (V and  components  u,  v),  energy  (e),  and  species  mass  fraction  (c,). 
Pressure  and  temperature  are  p  and  T,  whereas  L  is  a  characteristic  lengthscale  for  the  flowfield. 
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The  diffusion  coefficient,  D,  is  understood  to  be  the  multicomponent  coefficient,  Dim  (specie  i 
into  the  mixture),  which  is  related  to  the  binary  diffusion  coefficient,  Dj  (specie  i  diffusing  into 
specie/), 
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(2) 


where  Xt  is  the  mole  fraction  and  M(-  is  the  molecular  weight  of  species  i,  klD  is  the  collision 
integral  and  o,y  is  the  collision  diameter  (available  in  the  literature).  The  thermal  diffusion  ratio, 
kT,  is  assumed  to  be  1.0. 


The  relations  for  species  enthalpy,  mixture  energy  per  unit  volume,  and  pressure  are  given  by  the 
following  state  equations.  It  should  be  noted  here  that  appropriate  state  equations  for  a 
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high-temperature  plasma,  which  are  somewhat  more  complex  than  these,  are  at  the  same  time  not 
widely  available  in  the  literature.  The  present  choice  of  these  relations  will  affect  the  results. 


h  =  ei+Ri  T  ,  e  =  +  P{l'  +  V  )  +  V  p  cAH°  ,  p  =  p  RT  , 

r~  1  2  fi  ■>  y  y  La  M 


(3) 


where  A H°  is  the  fonnation  heat  for  species  i.  The  ratio  of  specific  heats  is  /  =  cp  /cy  ,  where  cp 
and  cv  are  the  mixture  specific  heats.  The  species  values  are  computed  using  polynomials  with 
coefficients  OC1-OC7  (22).  The  specific  gas  constant  is  R,  (the  universal  constant  is  R),  and  the 
heats  can  be  related  using  Cvi  =  Cpj-Ri . 


and 


N  N  f  Or  ^ 

cp  =  iicpi  ct  +  f  ^ 


i  =  l 


cpj  —  R:  +  a2T  +  +  a^T  +  a5T  +  +  cijT4  j 


(4) 


Likewise,  viscosity  and  thermal  conductivity  coefficients  are  given  by  polynomials  with 
coefficients  A,  B,  C,  and  D  given  in  reference  (22)  (coefficients  for  each  specie  are  different  for 
each  variable  and  each  specie).  The  mixture  values  of  viscosity  and  thermal  conductivity  are 
computed  using  empirical  mixture  rules  discussed  in  reference  (19). 

Mfi)  =  AlniT)  +  jr  +  jk  +  Di  ’  and  Mk)  =  A  HT)  +  jr  +  jk  +  A  •  (5) 


Chemical  reactions  can  be  expressed  in  a  general  reaction  equation.  X,  represents  the  symbol  for 
specie  i,  and  v  is  the  stoichiometric  coefficient  for  species  i  (primed  for  reactant  and 
double-primed  for  product)  in  the  reactions  (10).  The  chemical  source  terms  that  appear  in  the 
H  array  for  equation  1  (i.e.,  w)  are  computed  using  a  general  reaction  rate  equation.  The  reaction 
rates  for  each  reaction,  k  (“f  ’  for  forward  and  “b”  for  backward),  are  usually  expressed  in 
Arrhenius  fonn  (see  reference  [ 19]  for  details)  but  often  include  corrections  for  high-pressure 
effects  (see  reference  [70]  for  reaction  rate  data).  The  backward  rate  coefficients  are  computed 
from  the  forward  values  and  the  equilibrium  constant  for  each  reaction  (19). 


2>,X  i>;x.  and  w,  =Ml(v'l--v)\k,X{cf-kiY\Cl 

1=1  1=1  v  i  i 


(6) 


The  Navier-Stokes  equations  (equation  1)  are  written  in  integral  form  and  then  re-expressed  in  a 
semi-discrete  fashion  using  a  finite-volume  discretization  technique.  The  numerical 
computations  proceed  by  solving  the  semi-discrete  equation  on  each  computational  cell  using 
central  and  upwind  numerical  differences  along  with  flux-limiting.  Once  properly  discretized, 
the  resulting  set  of  algebraic  equations  is  solved  in  a  coupled  manner  in  time,  using  an  explicit 
time-accurate  method.  The  numerical  time-step  is  computed  using  the  CFL  condition.  A 
separate  chemical  time-step  is  computed  as  well.  The  final  time-step  is  based  on  the  smallest  of 
these.  Nusca  (19-21)  offers  a  complete  treatment  of  the  numerical  scheme. 
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5.  Results  and  Discussion 


The  computational  domain  chosen  to  be  used  to  simulate  the  experiment  and  sketched  in  figure  2 
extends  from  the  protective  plate  to  a  small  distance  along  the  pressure  probe  (0.1 8-m  total 
length)  and  from  the  centerline  of  the  capillary  (and  probe)  to  a  small  radial  distance 
(0.08-m  total  radius).  This  region  was  discretized  using  277  axial  grid  cells  and  160  radial  grid 
cells,  distributed  in  essentially  an  even  spacing  throughout.  The  boundary  conditions  for  this 
region  are  symmetry  on  the  axis  (Y  =  0),  outflow  at  X  =  0.18  m  and  Y  =  0.08  m, 
no-slip/no-penetration  on  the  plate  surface  (X  =  0,  0.0159  <  Y  <  0.08  m)  and  the  probe  surface, 
and  specified  inflow  at  the  capillary  exit  (X  =  0,  0  <  Y  <  0.0159  m).  The  inflow  conditions 
follow  from  those  displayed  in  figures  4-6  (note:  see  table  1  for  a  complete  list  of  species 
specified  at  the  inlet).  Initially,  the  entire  flowfield  is  filled  with  air  (0.8  mole  fraction  of  N2  and 
0.2  mole  fraction  of  O2). 

To  aid  in  the  interpretation  of  computational  results  for  the  plasma  jet,  figure  7  shows  a 
schematic  of  the  expected  gas  dynamic  features  as  have  been  observed  in  highly  underexpanded 
supersonic  jets  (25,  24).  The  efflux  of  plasma  from  the  inlet  generates  a  weak  precursor  shock 
(A)  that  expands  spherically.  Behind  this  shock  is  air;  the  plasma  is  entirely  contained  by  this 
shock  and  is  separated  from  the  air  by  an  irregularly  shaped  contact  surface  (B)  across  which 
pressure  and  velocity  are  preserved,  but  entropy  changes  discontinuously.  Expansion  waves  (C), 
generated  at  the  inlet,  travel  to  the  precursor  shock  (A),  are  reflected  as  weak  compression 
waves,  and  coalesce  into  a  strong  oblique  shock  or  barrel  shock  (D)  within  the  plasma  jet.  This 
barrel  shock  (D)  terminates  in  an  irregular  reflection  that  forms  a  triple-point  (E),  joining  the 
barrel  shock  (D),  it’s  reflection  (F),  and  a  normal  shock  (G)  or  Mach  disk.  Whereas  the 
precursor  shock  (A)  is  relatively  weak  and  diffuse,  producing  a  mildly  supersonic  flow,  the 
barrel  shock  (D)  and  Mach  disk  (G)  are  strong  shocks  that  enclose  a  fully  supersonic  flow 
region. 


Figure  7.  Schematic  of  gas  dynamic  features  in  a  highly  underexpanded  jet. 
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Figure  8  shows  the  time  history  of  the  plasma  jet  from  shortly  after  entrance  into  open-air 
(0.05  ms)  to  shortly  after  impact  upon  the  pressure  probe  (0.24  ms).  Note  that  the  frames  in 
figure  8  are  not  evenly  spaced  in  time  but  are  chosen  to  highlight  the  features  of  figure  7.  Mach 
number  contours  are  displayed  in  gray-scale  from  0  (bright  white)  to  1.6  and  above  (dark  black). 
As  the  low-density  (0.08-kg/m  )  plasma  efflux  enters  the  open-air,  a  precursor  shock  is  formed. 
At  this  time,  the  jet  is  moving  at  -5600  m/s,  which  is  -Mach  1.6,  based  on  the  local  sound  speed 
in  the  mixture  (note  that  the  plasma  gas  is  of  very  low  molecular  weight,  -4.3).  By  -0.20  ms, 
the  predominant  gas  dynamic  features  in  the  jet  have  been  established,  including  the  intersection 
of  the  barrel  shock  with  the  Mach  disk  (X  =  0.06  m)  and  a  weak  reflected  shock.  Due  to  the 
variable  viscosity  in  the  flowfield  (see  equation  5),  the  precursor  shock  is  more  diffuse  than  the 
nonnal  shock  of  the  Mach  disk  or  the  oblique  “barrel”  shock.  Note  the  subsonic  flow  region 
between  the  normal  shock  and  the  precursor  shock;  within  this  region  resides  the  contact  surface. 
The  precursor  shock  reaches  the  pressure  probe  at  -0.22  ms.  Stagnation  of  supersonic  flow  on 
the  probe  tip  causes  a  normal  shock  to  standoff  from  the  probe.  By  0.24  ms,  this  tip  shock  is  no 
longer  required  because  the  flow  is  subsonic  in  the  probe’s  near  field.  Figure  9  shows  the 
measured  and  computed  pressure  at  the  probe  tip.  As  was  discussed  in  relation  to  figure  3,  the 
computed  pressure  history  has  been  shifted  by  0.027  ms.  Even  so,  the  computed  precursor  shock 
arrives  at  the  probe  -0.02  ms  early  and  registers  a  pressure  that  is  -30%  higher  than  measured. 
Computed  pressure  oscillations  around  the  peak  are  due  to  the  unsteady  nature  of  the  normal 
shock  at  the  probe  tip.  The  presence  of  particulates  in  the  flow  stream  reported  by 
White  et  al.  (5)  may  account  for  some  of  this  discrepancy  as  may  the  absence  of  both  a  true 
plasma  equation  of  state  (see  equation  3)  and  radiation  energy  losses  in  the  computations.  At 
this  point  in  time,  the  computation  was  halted  until  further  investigations  of  the  source  of  the 
pressure  overprediction  was  discovered. 

Figures  10  and  1 1  show  computed  results  for  the  plasma  jet  at  0.13  and  0.22  ms,  respectively,  in 
terms  of  eight  variables.  The  pressure  and  temperature  have  been  nonnalized  by  one  atmosphere 
and  300  K,  the  values  ahead  of  the  precursor  shock,  respectively.  The  heat  capacity  ( cp ),  thermal 
conductivity  (k),  and  viscosity  (//)  have  been  nondimensionalized  using  the  values  for  air; 
therefore,  ahead  of  the  precursor  shock,  these  variables  are  1.0.  At  0.13  ms  (figure  10),  the 
precursor  shock  is  indicated  by  the  sharp  pressure  gradient  at  X  =  0.09  m  from  the  inlet  (the 
spherical  shape  of  this  shock  is  distorted  by  the  plotting  scale).  The  nonnal  shock  (Mach  disk)  at 
X  =  0.035  m  is  clearly  visible  in  the  pressure  contour  map,  as  is  the  high-pressure  flow  from  the 
inlet.  Behind  the  precursor  shock  and  in  front  of  the  nonnal  shock  (i.e.,  0.035  <  X  <  0.09  m),  the 
pressure  and  temperature  are  quite  high,  indicating  extreme  compression  of  the  flow  due  to  these 
two  shocks.  The  plots  of  heat  capacity  and  thermal  conductivity  highlight  the  region  of  plasma 
(note  the  piston-like  shape  with  expanded  leading  edge)  in  which  these  variables  take  on  values 
that  are  much  higher  than  those  in  ambient  air.  The  plot  of  viscosity  indicates  that  the  plasma 
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Figure  9.  Measured  and  computed  pressure  at  the  probe  tip. 

itself  represents  largely  laminar  flow  (i.e.,  high  value  of  viscosity).  The  plot  of  hydrogen  (H) 
atom  concentration  confirms  that  the  plasma  consists  mainly  of  H  gas  (figure  6)  and  electrons 
(the  plot  of  electron  concentration  is  qualitatively  similar  to  the  plot  of  H).  In  the  high-pressure 
and  high-temperature  region,  the  density  is  also  high  and  causes  the  largest  local  concentration  of 
H  away  from  the  inlet.  The  formation  of  OH  and  CO?  will  not  take  place  until  sufficient  thennal 
and  chemical  conditions  arise  in  the  flowfield,  as  dictated  by  the  detailed  chemical  kinetics 
mechanism  (the  reader  is  referred  to  reference  [1 0]).  The  contour  plots  of  these  species  indicate 
that  such  conditions  do  exist  outside  the  plasma  core  region  (0  <  X  <  0.055  m,  0  <  Y  <  0.018  m 
at  this  time)  where  mixing  with  air  (i.e.,  O2  and  N2)  is  sufficient.  Space  limitations  do  not  permit 
the  displaying  of  all  chemical  species  distributions,  but  this  information  is  being  probed  as  a 
means  of  characterizing  the  chemical  nature  of  the  plasma  efflux.  Figure  1 1  shows  the  same 
variables  plotted  for  0.22  ms  in  time,  when  the  leading  edge  of  the  plasma  jet  (i.e.,  the  precursor) 
has  just  reached  the  probe.  Flow  stagnation  on  the  probe  tip  causes  local  compression,  and  a 
bow  shock  is  generated  on  the  tip  (see  pressure  plot).  Gas  temperature  (and  density)  also 
increases  near  the  probe  as  does  the  concentration  of  H.  Flow  heating,  species  mixing  with  air, 
and  chemical  reactions  cause  increased  concentrations  of  OH  (especially)  and  CO2;  this  indicates 
increased  levels  of  chemical  reaction  in  the  flowfield.  Of  course,  the  probe’s  near-field 
properties  continually  change  as  subsequent  regions  of  the  jet  flow  impinge  upon  the  tip 
(i.e.,  after  the  precursor  shock,  the  contact  surface,  and  the  Mach  disk  containing  a  region  of 
supersonic  flow). 
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Figure  10.  Computational  results  displayed  in  eight  contour  variables  for  0.13  ms. 
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After  looking  at  these  results,  one  can  suppose  that  the  physical  and  chemical  nature  of  the 
plasma  efflux  in  the  vicinity  of  the  probe  is  very  different  from  that  near  the  capillary.  Gas 
dynamic  expansion  in  open  air  has  caused  the  generation  of  various  shocks  and  mixing  layers 
that  have  combined  the  plasma  species  with  air  and  have  resulted  in  chemical  reactions.  Table  1 
lists  the  species  considered  in  the  present  model  and  provides  the  concentrations  for  a  time  of 
0.22  ms  (i.e.,  displayed  in  figure  1 1)  and  at  three  locations  along  the  centerline  of  the  jet:  at  the 
inlet  (i.e.,  capillary  exit  plane),  just  outside  of  the  inlet,  just  ahead  of  the  Mach  disk,  and  within 
the  precursor  shock  at  the  stagnation  region  on  the  probe  tip.  It  is  interesting  to  note  the 
increased  concentrations  of  various  species  at  the  probe  tip,  especially  those  species  that  are  not 
predominant  at  the  capillary  (CO,  NO+,  N2,  O,  0+,  OH,  H20,  C02),  and  those  abundant  in  the 
plasma  that  have  collected  on  the  tip  (H,  H+,  C+,  electrons).  The  C  atom,  abundant  in  the  plasma 
efflux  at  the  capillary,  is  absent  at  the  probe  tip  (possibly  combined  in  CO  and  C02).  These 
results  confirm  the  supposition.  The  ultimate  goal  is  to  understand  how  the  propellant  responds 
to  the  plasma  jet  both  mechanically  (i.e.,  pressure  and  thermal  loading  from  the  shock  system) 
and  chemically.  This  information  concerning  the  pressure,  temperature,  and  chemical  conditions 
at  the  probe  tip  are  necessary  starting  conditions  for  subsequent  modeling  of  solid  propellant 
combustion  in  the  plasma  stream. 


6.  Conclusions 


A  time-accurate  CFD  code  has  been  applied  to  the  modeling  of  the  high-temperature,  chemically 
reactive  plasma  efflux  from  an  ETC  igniter  fired  into  open  air.  The  CFD  code  has  been  linked  to 
a  plasma  capillary  model  for  purposes  of  specification  of  the  exit  conditions  of  the  igniter.  The 
major  features  of  this  efflux  have  been  resolved  by  numerical  simulation  revealing  a  highly 
underexpanded  jet  with  a  strong  precursor  shock,  a  barrel  shock  that  reflects  at  a  triple-point,  and 
a  normal  shock  or  Mach  disk.  Impact  of  the  jet  upon  a  probe  generates  an  intermittent  bow 
shock  on  the  probe  tip  that  is  dissipated  in  the  low-velocity  flow  ahead  of  the  Mach  disk  and  is 
then  re-established  with  the  passage  of  this  nonnal  shock.  The  presence  of  this  shock  system  has 
important  implications  for  mechanical  damage  to  solid  propellant  exposed  to  the  plasma  efflux. 
Chemical  conditions  at  the  probe  are  quite  different  from  those  at  the  ETC  igniter.  These 
computer  simulations  represent  the  first  such  detailed  modeling  of  ETC  igniter  plasma  efflux. 
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